Direct conversion of rheological compliance measurements into storage and loss moduli 
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We remove the need for Laplace/inverse-Laplace transformations of experimental data, by pre- 
senting a direct and straightforward mathematical procedure for obtaining frequency-dependent 
storage and loss moduli (G'{u)) and G"{ijj) respectively), from time-dependent experimental mea- 
surements. The procedure is applicable to ordinary rheological creep (stress-step) measurements, as 
well as all microrheological techniques, whether they access a Brownian mean-square displacement, 
or a forced compliance. Data can be substituted directly into our simple formula, thus eliminating 
traditional fitting and smoothing procedures that disguise relevant experimental noise. 



The procedure that has become established [H , for ob- 
taining frequency-dependent dynamic moduh from non- 
oscillatory rheometry, is to fit the experimental data to a 
particular model (often the generalized Maxwell model is 
used), and subsequently to calculate the resulting com- 
plex viscoelastic modulus, G*{oj) = G'{lj) + iG"{uj), for 
that parametrization of the model. That procedure can 
be somewhat restrictive, as it may force the user to ap- 
proximate their data into the prescribed form, or to use 
a very large number of fitting parameters. It also artifi- 
cially hides experimental noise, making the uncertainties 
in the final results difficult to quantify. Equivalently, one 
can find an approximate Laplace transform of the time- 
dependent data , then derive the Laplace transform of 
the stress relaxation modulus, and subsequently trans- 
form from a Laplace to a Fourier description (either nu- 
merically or, for certain functional forms, analytically). 
In either case, the procedure limits the user's freedom in 
the types of formulae that can easily be fitted and manip- 
ulated, and it can be somewhat laborious. Furthermore, 
in cases where the experimental data are imperfectly fit- 
ted by the pre-conceived functions, or other approxima- 
tions are introduced (such as an approximate Laplace 
transform or its inverse Q), the accuracy of the de- 
rived moduli becomes vague. Here, we show that a more 
direct, straightforward and accurate treatment of rheo- 
metric data is possible, and derive a formula for G*{lu) 
in terms of the experimental data points themselves. 

Let us begin by summarising why the conversion of 
time-dependent rheometry data into storage and loss 
moduli is traditionally such a complicated process. (For 
more discussion and approximate solutions, see 0, [H].) 
In principle, a simple relationship exists between the dy- 
namic moduh, G'{uj) and G"{ui), and the experimentally 
accessible time-dependent compliance, J{t). In a shear- 
creep (or stress-step) experiment, 



J{t) = 7(i)/ao 



(1) 



is the ratio of the time-dependent shear strain jit) to the 
magnitude ctq of the constant stress that is switched on 
at time t = 0. Alternatively, in a passive microrheology 
experiment, the function J(t) is proportional to the mean 
square displacement of a probe particle 0, executing 
Brownian motion in the viscoelastic fluid (averaged over 
many trajectories). Since the relaxation modulus is re- 
lated to the compliance by a convolution , 



G{t) J{t - t) dt 



(2) 



it is, in principle, a simple matter to extract the modu- 
lus by deconvolving Eq.Q, using an integral transform, 
such as the Fourier transform. So the Complex viscoelas- 
tic modulus G*{lj) (which is the Fourier transform of 
the time derivative of G{t)), is a simple function of the 
Fourier-transformed compliance: 



G*iuj) 



1 



iujj{uj) 



(3) 



where iujj = J* is sometimes called the dynamic compli- 
ance The problem that arises is due to the fact that 
the Fourier transform (denoted !F[. 
ance. 



.](w) ) of the compli- 



J{u) ^ T[Jluj) 



J{t) e-*"* dt. 



(4) 
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is not a convergent integral, since J(t) grows with increas- 
ing time. (Even for a solid, where J(t) tends to a finite 
constant at long times, the integral in Eq. Q remains 
undefined.) This has led investigators to resort instead 
to a Laplace transform, Jl{s) = J(i) e~** dt, since 
its integral is convergent. Deconvolution of Eq. ^ thus 
yields Gl{s), the Laplace tranform of G{t). Only if the 
resulting function is expressed as a simple formula, then 
the desired complex modulus can finally be obtained by 
analytic continuation, G*{uj) — s Gl{s)\^^,-^. 

In this Letter, we discuss how to bypass the above foray 
into Laplace space, by working directly with the Fourier 
transform of the compliance. Although its integral repre- 
sentation in Eq. ^ is not convergent, the quantity J{uj) 



is nonetheless well defined, as is apparent in Eq. Q , since 
G*{ijj) exists for all real finite to. It is well known how 
to find the Fourier transform of an unbounded function 
such as J{t). We shall nevertheless introduce the method 
in pedagogical detail, in order to clarify its applicability 
to any causality-respecting compliance function, includ- 
ing that dictated by the raw experimental data. This 
will allow us to find a simple formula for the viscoelastic 
moduli that are implied by the data. 
First we note that causality requires 



J{t) =0 for t<0 



(5) 



as there can be no response before the stress step is ap- 
plied. So J{t) is a function resembling the sketch in 
Fig. [T^, and J{uj), required in Eq. ([3]), is its Fourier 
transform. In the long-time limit, the response of a fluid 
(viscoelastic or otherwise) to an imposed step stress is to 
undergo shear at a constant rate. So that the compliance 
J(t), which is proportional to the total strain, asymptotes 
to a straight line (see Fig. [TJi) with a gradient equal to 
the reciprocal of the static viscosity. So J{t), the second 
derivative of J(t), is a function that vanishes at large t, 
and its Fourier transform therefore converges. We can 
reconstruct the former Fourier transform from the latter, 
since they are simply related thus: 



Jiu;) 



-1 



(6) 



It is not immediately obvious that a function can be re- 
constructed from its second derivative, in this way, as 
there is a danger of losing information about the absolute 
offset and slope of the original function J{t). However, we 
can retain that information since we have, as a reference, 
the known part of the function given in Eq. ([5|). To do so, 
we must perform the double differentiation of J{t) over 
the whole of the function's domain, including negative 
and zero values of t, as shown schematically in Fig. [T]3 
and c. Notice that J{t) has a discontinuity at t = 0, of 
size J(0), the initial gradient (at time t = 0+) of the com- 
pliance. Differentiating the discontinuous function J(t) 
yields the second derivative, J{t) which vanishes at neg- 
ative times, has a Dirac delta function of strength J(0) 
located a.t t — 0, and is finite for positive t (see Fig. [TJ;). 
Substituting that form into Eq. ^ yields 



-j(0) 



1 

./o+ 



■ J{t) dt, 



(7) 



where the integration is over positive values of t only, for 
which the function J(t) is non-singular. That integration 
can be performed by standard Fourier transform methods 
if we define a simpler function J2{t) that is the second 
derivative of J{t) excluding the delta function. 



Mt) 



Jit) 



for t > 0, 
otherwise. 



(8) 




(c) 
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FIG. 1: (a) Sketch of a typical time-dependent compliance, 
J{t), which must vanish for negative t due to causality, (b) 
Its first derivative, J{t). (c) The second derivative, J{t). 



In terms of J2(t), the integral in Eq. ([7]) is simply 
J-[J2]{oj), so that Eq. ([3]) yields a simple formula for the 
storage and loss moduli, G'{uj) + iG"{uj) — iuj/{J{0) + 

Let us generalize to include compliance functions that 
are discontinuous at t = 0, since a finite discontinu- 
ity J(0) = limt^o+ J(t) 7^ is often observed, reflect- 
ing the fact that the data-acquisition rate cannot ac- 
cess the regime of a material's response preceeding a 
small-f plateau in J. Such a compliance function has 
a delta- function contribution to its first derivative J{t), 
of strength J(0). Hence, defining Ji{t) to exclude that 
delta function, so that J{t) = Ji{t) + J{0)6{t), we have 
J{uj) = {J(0) + T[Ji]{u;)} / iiv , ultimately yielding 



G'{Lu) + iG"iuj) = 



iujJ{Q) + J{0)+T[J2]{uj) 



(9) 



Since Eq. ([9]) holds for any compliance function, we 
now apply it directly to the experimental data, by defin- 
ing the piecewise linear function J{t) that interpolates 
between data points, depicted in Fig. [2l It is obvious 
that the resulting angularity and non-monotonicity of 
the function are unphysical consequences of the exper- 
imental noise, but attempting to remove them by using 
a smooth fitting function would constitute doctoring the 
data. Instead, we take the data at face value by us- 
ing this function that passes through every data point. 
Since evaluation of the Fourier transform requires J{t) 
to be defined for all positive t, while experimental data 
are finite, our piecewise linear function must include an 
extrapolation to t = oo, thus introducing one extra pa- 
rameter, the steady-state viscosity 77 (see Fig. [2]). This is 
not a peculiarity of the present method; any data-analysis 
for converting J{t) to G*{uj) requires such extrapolation, 
though some methods obscure it in arcane algorithms. 
The second derivative of our experimental function J{t) 
is a series of delta functions, so its Fourier transform is 
trivial to evaluate. The strengths of the delta functions 
are equal to the discontinuities in gradient of J(t), and 
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J(0) 



FIG. 2: An example of data points {tk, Jk) from experimental 
measurements of the compliance as a function of time. The 
data are interpolated by a piecewise-linear function, and also 
linearly extrapolated to infinity by a final line segment of 
gradient »7~^, which vanishes in the case of elastic solids for 
which the present treatment is equally valid. 



Eq. ([9]) becomes straighforwardly expressed in terms of 
the experimental data points {ti, Ji), which need not be 
equally spaced: 

^ = ..j(o) + (i-.-'.)i^« + ::^ 

G*{uj) w V ' h rj 



N 



k=2 



Jk — Jk- 



tk — tk~l 



j2 ( -^^^ ) (e— _ e— (10) 



Note that Eq. pH)) requires the convention > so 
that, if the compliance is non-zero at t = 0, it enters the 
formula in the value of J(0) only. 

Equation (jlOp makes a direct hnk between the exper- 
iment and the resulting graphs of G" and G", thus re- 
moving any subjective judgement from the results, and 
allowing genuine experimental noise and uncertainties to 
appear on those graphs for critical evaluation. The for- 
mula also has the advantage of being very easy to ap- 
ply. For instance, the code required to evaluate it us- 
ing Mathematica® can be written in just three lines, as 
follows i- 

{t , J}=Transpose [Import ["filename .txt" , "Table"] ] ; 

G[lu_, J0_, ri_, {t_, J_}] := I uj/(I LU 30 + 
(1-Exp[-I LJ t[[l]]])(J[[l]]-JO)/t[[l]] + 
Exp[-I LU t [ [Length [t]]] ] / + 
Sum[ (Exp[-I LU t[[k-l]]]-Exp[-I uu t [ [k] ] ] ) 
(J[[k]]-J[[k-l]])/(t[[k]]-t[[k-l]]), 
{k, 2, Length [t]} ] ) ; 

LogLogPlot [{Re [Gluj,0. 0000023 , 1 145300 , {t , J}] ] , 

Im [G [cj , . 0000023 ,1145300 , {t , J}] ] } , , . 001 , 1000}] 

(Some installations of Mathematica® first require a li- 
brary to be loaded, using << Graphics ' Graphics ' .) 
The first line imports the experimental data (a 
list of pairs of numbers) from a file (here named 
f ilensune . txt). The second line defines the function 
in Eq. (fTO|) . and the third displays the resulting real 
and imaginary parts of the complex modulus, in this 
case using the parameter values J(0) = 2.3 x 10~®, 
rj = 1.1453 X 10® that characterise our data. 
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FIG. 3: Creep comphance J vs. time t, for PI (M„ = 150000, 
M^/Mn = 1.03) measured at T = 0°C. Top left inset: 
strain vs. time for three different applied stresses of magni- 
tude 200 Pa (green triangles), 500 Pa (red circles) and 1000 Pa 
(black squares). Bottom right inset: compliance vs. time, 
for the same data sets. The data-collapse testifies to the lin- 
earity of the creep measurements. 



We have used Eq. (|T0|) to obtain the frequency- 
dependent storage and loss moduli of a near- 
monodisperse polyisoprene melt with a weight-average 
molar mass Mw of 152 kg/mol and polydispersity 
Mw/Mn of 1.03 (where Mn is number-average mo- 
lar mass). The time-dependent creep compliance J(t) 
was determined in a stress-step measurement performed 
at Celsius using a commercial AntonPaar MCR-501 
rheometer (cone diameter 25mm, cone angle 1 degree). 
The insets of Fig. [3] show the measured strain curves 
(top left inset), and thus the compliance curves (bottom 
right inset) , at three different shear stresses of magnitude 
200 Pa, 500 Pa and 1000 Pa. The good superposing of 
the three compliance curves indicates that the applied 
stresses were small enough to access the fluid's linear 
regime (where the ratio on the RHS of Eq. ^ is inde- 
pendent of (To) J but large enough to provide a satisfactory 
signal-to-noise ratio. 

The data points in Fig. [3] (main graph) are the average 
of three compliance measurements. Scatter due to ex- 
perimental noise is apparent, particularly at small t, but 
these raw data were substituted directly into Eq. (fTO]) 
without smoothing or fitting. The resulting functions 
G'{lu) and G" [u) are plotted in Fig. [5] (red and blue 
curves respectively). Notice that the curves are not 
smooth. This demonstrates a virtue of our straightfor- 
ward deconvolution method (Eq. (UU])) over established 
methods: that it preserves the experimental noise. The 
noise visible in the curves is a true refiection of the ex- 
perimental uncertainties. 

For comparison, Fig.|4]also shows data from oscillatory 
measurements on the same fluid. The agreement is good 
across flve orders of magnitude in modulus and in fre- 
quency. Not only are features such as the characteristic 
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FIG. 4: Storage (red solid line) and loss (blue dotted line) 
moduli, respectively the real and imaginary parts of G* 
found by substituting compliance data (Fig. [S} into Eq. (|10|l . 
Crosses are data from oscillatory rheometry of the same fluid, 
shown for comparison. Inset: To demonstrate the efficacy 
of Eq. (|10|l for smooth functions also, the compliance data 
(Fig. O were fitted to the function J{t) — (t/ri) + a + 
b tanh (c + d In t + e In^ t + f In'^ , the form of which is mo- 
tivated only by fit quality. The seven parameter values were 
77 = 1.1453 X lO'^Pa, a = 4.7 X 10"'^Pa-\ 6 = 2.4 x 10"^Pa-\ 
c = 0.288082, d = 0.247501, e = 0.0174205, / = 0.000685812, 
yielding J(0) = 2.3 x lO'^'Pa"^ The result of substituting 
200 sampled points from this function, at equal logarithmic- 
time intervals, into Eq. (|10|l are plotted in the inset, with the 
oscillatory data again reproduced for comparison. 

relaxation times (e.g. where the moduli cross) accurately 
obtained, but absolute values of the moduli are consis- 
tent, up to random errors in experimental reproducibil- 
ity, thus demonstrating that the same linear rheological 
information can be extracted from the stress-step exper- 
iment as from a large number of oscillatory experiments. 

The moduli from Eq. (fTO)) have been plotted in Fig. [4] 
for angular frequencies uj in the domain uj^i^ < uj < 
Wmax- Outside this frequency window, the moduli given 
by Eq. pO|) are dominated by artifacts. The lowest ac- 
cessible frequency, Wmin ~ ^/v^' determined by the ex- 



periment's duration t^r. Obtaining a single data point by 
oscillatory rheometry at that same frequency Wmin would 
require several complete oscillations, thus taking an or- 
der of magnitude longer than the creep measurement of 
the entire dynamic spectrum. The highest accessible fre- 
quency, cjmax ~ ^r^' early-time resolution of 
the stress-step experiment, where the first reliable data 
are obtained at time ti. At high frequencies, however, 
oscillatory measurements are relatively quick to perform 
and can yield better precision than can be obtained from 
the short-time creep response of the rheometer (depend- 
ing on details of the instrument's design). Hence, in prac- 
tice, using a combination of creep and oscillatory mea- 
surements may be the best strategy to determine a fluid's 
entire dynamic spectrum. 

Finally, we note that Eq. PH)) has another use be- 
sides substitution of raw data. Even for analytical func- 
tions J(t) (such as might be used in theoretical work, 
or to approximate noisy data), standard numerical al- 
gorithms can fail to evaluate the Fourier transform re- 
quired in Eq. ([5]), if the small-t behaviour of J is non- 
trivial. Equation (|10p turns out to be a reliable method 
for numerically evaluating the required function, with 
accuracy greatly superior to simple quadrature algo- 
rithms such as trapesium rule. To demonstrate this, a 
smooth function (see caption to Fig. [3]) was fitted to 
the compliance data in Fig. [31 and was then sampled 
at logarithmically-uniform intervals in t, for substitution 
into Eq. (|10p . This yielded the smooth curves in the in- 
set to Fig. m We found convergence on the exact result 
for 200 sample points, and that the above Mathematica® 
code executed in a matter of seconds, whereas the built- 
in FourierTransf orm algorithm stalled when applied to 
evaluate Eq. 

In summary, dynamic moduli can be straightfor- 
wardly obtained by substitution of compliance data into 
Eq. (Uni), which is equally valid for viscoelastic fluids or 
solids. The equation is quick to evaluate, removes the 
need for approximate fitting or obscure black-box algo- 
rithms, and correctly preserves the experimental noise 
that is so crucial to good scientific methodology. 
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